Slow cross-symmetry phase relaxation in complex collisions 

Luis Benet, 1 Lewis T. Chadderton, 2 Sergey Yu. Kun, 3 ' 4 '[*| Oleg K. Vorov, 5 and Wang Qi 1 

1 Instituto de Ciencias Fisicas, Universidad Nacional Autdnoma 
de Mexico (UNAM), 62210-Cuernavaca (Morelos), Mexico 
2 Atomic and Molecular Physics Laboratary, RSPhysSE, 
Australian National University, Canberra ACT 0200, Australia 
3 Facultad de Ciencias, Universidad Autdnoma del 
Estado de Morelos, 62209-Cuernavaca (Morelos), Mexico 
^Nonlinear Physics Center and Department of Theoretical Physics, 
RSPhysSE, Australian National University, Canberra ACT 0200, Australia 
5 Department of Physics and Astronomy, 
Drake University, Des Moines, Iowa 50311, USA 
6 Institute of Modern Physics, Chinese Academy of Sciences, Lanzhou 730000, China 

We discuss the effect of slow phase relaxation and the spin off-diagonal S-matrix 
correlations on the cross section energy oscillations and the time evolution of the 
highly excited intermediate systems formed in complex collisions. Such deformed 
intermediate complexes with strongly overlapping resonances can be formed in heavy 
ion collisions, bimolecular chemical reactions and atomic cluster collisions. The 
effects of quasiperiodic energy dependence of the cross sections, coherent rotation 
of the hyperdeformed ~ (3 : 1) intermediate complex, Schrodinger cat states and 
quantum-classical transition are studied for 24 Mg+ 28 Si heavy ion scattering. 

PACS numbers: 25. 70.-z, 24.10. Cn,03.65.Nk,03.65.-w 



I. INTRODUCTION 

The dominating idea in the modern theory of highly excited strongly interacting systems 
is that phase randomization time is the shortest time scale of the problem Applica- 
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tion of this idea to the theory of quantum chaotic scattering for colliding systems with 
rotationally invariant Hamiltonians implies the absence of correlations between the reaction 
amplitudes carrying different total spins lj. However, though seemingly plausible, the as- 
sumption of a very fast phase relaxation is not consistent with many data sets on complex 
quantum collisions. In particular, the anomalously long-lived spin off-diagonal ^-matrix 
correlations have been identified from the data on forward peaking of evaporating protons 
in nucleon induced 

as and photonuclear [4| reactions. Such long-lived correlations reflect 
an anomalously slow phase relaxation, which is many orders of magnitude longer than the 
energy relaxation. This provides a manifestation of a new form of matter: thermalized 

n n 

non-equilibrated matter introduced by one of us in Refs. [5|, 16|]. The effect is of primary 
importance for many-qubit quantum computation since anomalously long "phase memory" 
can extend the time for quantum computing far beyond the quantum chaos border 0, Q] . 
An effect of a very slow phase relaxation has also been strongly supported by numerical 
calculations for H+D2 [7], F+HD [8] and He+H^ state-to-state chemical reactions. In 
these calculations, a slow phase relaxation manifests itself in stable rotating wave packets 
of the intermediate complexes [lO]. Interestingly, this same effect of stable coherent rota- 
tion was originally revealed for heavy ion collisions, e.g., for 19 F+ 89 Y ll|], 28 Si+ 64 Ni 12], 



12/ 



C+ 24 Mg [lahil, 24 Mg+ 24 Mg and 28 Si+ 28 Si [15|, 58 Ni+ 46 Ti and 58 Ni+ 62 Ni [To) collisions. 



In this paper we reveal the effect of slow phase relaxation for yet another heavy ion 
scattering system, 24 Mg+ 28 Si. This effect is studied in relation to quasiperiodic energy de- 
pendence of the cross sections, coherent rotation of the hyperdeformed ~ (3 : 1) intermediate 
complex, Schrodinger cat states and quantum-classical transition for 24 Mg+ 28 Si heavy ion 
scattering. 



II. CROSS SECTION ENERGY AUTOCORRELATION FUNCTION 

We consider first the scattering of spinless collision partners with spinless collision fragments 
in the exit channel. Using the semiclassical asymptotics of Legendre polynomials for J 3> 1, 
where J is the total spin of the system, we represent the cross section in the form 



da(E, 9)1 d9 = a(E, 6) = a d (6) + 5a(E, 6), 



(1) 
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with 5a(E, 9) = 5aM{E,6) + 5<T^{E,6), Sa^(E,9) = \Sf^(E,9)\ 2 and 

5fW(E, 9) = ]T(2J + 1)W(J) 1/2 5S J (E) exp[iJ($ ± 9)}. (2) 
j 

In these expressions, <Td(9) = \F d {9)\ 2 is the energy independent potential scattering cross 
section, $ is the average deflection angle obtained from a linear approximation for the J- 
dependence of the potential phase shifts in the entrance and exit channels, W(J) is the 
average partial reaction probability, and 5S J (E) are normalized, (\5S J (E) | 2 ) = 1, energy 
fluctuating around zero ^-matrix elements corresponding to time-delayed collision processes. 
The brackets (...) stand for the energy E averaging. In the expression for a(E,9) we 
have dropped (i) the highly oscillating angle interference term between the 8f^ + \E,9) and 
5f(~~\E, 9) amplitudes, and (ii) the interference terms between the energy smooth potential 
scattering amplitude Fd(9) and energy fluctuating amplitudes 5f^(E,9). This is because 
the excitation functions data for 24 Mg+ 28 Si scattering 17J] were obtained by averaging over 
a wide A0 c . m . ~ 77° — 98° angular range. 

In calculating the cross section energy autocorrelation function, 

C(e) = (a(E + e, 9)a(E, 9))/(a(E, 9)) 2 - 1, (3) 



we take into account the 5-matrix spin off-diagonal correlation 181 ] 



(SS J {E + e)5S J (E)*) = T/{T + (3\J-J'\+ ihu{J - J') - is). (4) 

Here, u is the angular velocity of the coherent rotation of the intermediate complex, {3 is 
the spin phase relaxation width and V is the total decay width of the intermediate complex. 

We take W(J) in the J-window form, W(J) = W(\J — I(E)\/g), where the average 
spin 1(E), for a given cm. energy of the collision partners, is close to the grazing orbital 
momentum. The J-window width g relates to the effective number of partial waves, g + 1, 
contributing to Sa(E,9). For the analyzed 24 Mg+ 28 Si scattering we estimate g — 1 — 5 , 



which is revealed by the shape of the measured elastic scattering angular distributions [17] 
corresponding to the maxima of the excitation functions. Although these angular distri- 
butions show regular oscillations with a well-defined period, they clearly deviate from the 
square of a single Legendre polynomial. We estimate the energy dependence of 1(E) in the 
linear approximation and obtain 1(E) = I + I(E — E) / AE . Here, J = 1(E), E is the energy 
corresponding to the center of the energy interval over which the cross section is measured, 
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AE = 2(E — B)/I and B is the Coulomb barrier for the collision partners in the entrance 
channel for a configuration of the two touching spherical nuclei 24 Mg and 28 Si. 

We calculate C{e) under the conditions g > 1, (3 < T and (3 <C Two. We take W( J) in the 
Gaussian form, W(J) oc exp[— (J — 1(E)) 2 / g 2 }. Generalizing the calculations in 19( to the 
case of finite f3 and arbitrary AE for the normalized (C(e = 0) = 1) cross-section energy 
autocorrelation function, we obtain 

= exp[-e 2 /2(huj) 2 d 2 } r exp[i7r\e\/(Tw - i(3)] i 

[£) Re[l/[l-exp[-7rr/(^a;-i/3)]]] W - exj>[iir(\e\ + iT) / (Two - ifi)}^ 1 j 

where d 2 = g 2 /(l — hu/AE) 2 . The above expression for C(e) has been obtained for d > 1. 
One can see that (3/h has the physical meaning of the imaginary part of the angular velocity 
signifying the space-time derealization of the nuclear molecule and the damping of the 
coherent rotation of the intermediate complex 18(. For (3 = 0, C(e) is an oscillating periodic 
function with period Two. For finite (3, the amplitude of the oscillations in C(e) decreases 
with increasing |e|: the larger (3 the stronger the damping of the oscillations. For (3 = and 
hoo = AE, Eq. (5) transforms to the result in Ref. 19]. 

Although Eq. (5) is obtained for spinless reaction fragments it also holds for reaction 
products having intrinsic spins. This can be shown using the helicity representation for the 
scattering amplitude 181 ] . 

In Fig. 1 we present the energy autocorrelation functions for the 24 Mg+ 28 Si elastic and 
inelastic scattering [20], constructed from the data on the excitation functions measured on 
the E c . m . = 49 — 57 MeV energy interval [17|. The experimental excitation functions were 
averaged on the 8 cm = 77° — 98° angle interval. The experimental C(e)'s are not Lorentzian 
but oscillate with a period ~ 0.75 MeV. The fit of the experimental C(e)'s for all the elastic 
and inelastic channels is obtained with T = 0.15 MeV, (3 = 0.1 MeV, hu = 0.75 MeV and 
d = 5. The calculated C(e)'s are normalized to the experimental data at e = 0. 

The extracted value of hu suggests an anomalously strong deformation of the interme- 
diate complex. Indeed, for J ~ 34 — 38 [17j, using the moment of inertia of a ~ (2 : 1) 
superdeformed intermediate complex, corresponding to the two touched spherical colliding 
nuclei 24 Mg and 28 Si, we have hu ~ 1.9 MeV. This value is bigger by a factor of about 2.5 
than the period of oscillations in the experimental C(e)'s. This reveals the excitation of 
~ (3 : 1) hyperdeformed coherent rotational states of the intermediate complex. 

It should be noted that the intrinsic excitation energy of the intermediate complex, which 
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we obtain by substracting deformation and rotation energy for the total energy, is about 
15 MeV or more. This corresponds to the average level spacing of D ~ 10~ 6 MeV or less. 
Therefore, the intermediate complex is in the regime of strongly overlapping resonances, 
Y/D > 10 5 . In this regime, the theory of quantum chaotic scattering and random matrix 
theory are conventionally assumed to apply [l| . In accordance with these approaches, which 
in particular reconfirmed the Ericson theory of the compound nucleus cross-section fluctua- 



tions 2ll | . the spin off-diagonal ^-matrix correlations vanish yielding C(e) = 1/[1 + (e/Y) 2 ]. 
The Lorenzian curves presented in Fig. 1 with Y = 0.85 MeV to fit the experimental data 
at < 0.1 MeV are clearly in contrast with the oscillations in all the experimental C(e)'s. 
Within our approach, the limit of vanishing spin off-diagonal S'-matrix correlations corre- 
sponds to (3 3> r, where %/(3 is the characteristic spin phase relaxation time. Therefore, 
the persistence of the oscillations in C(e) indicates an anomalously long spin off-diagonal 
"phase memory" . 

In Fig. 1 we present another possible fit of the experimental C(e)'s with the same T = 0.15 
MeV and Two = 0.75 MeV, but with the different values (3 = 0.03 MeV and d — 1. One can 
see that both the fits are qualitatively and quantitatively undistinguishable. The question 
arises if the quantities of the interest, in particular the phase relaxation width /3, can reliably 
be determined from the data. 

III. TIME POWER SPECTRUM OF THE COLLISION 

Consider the time (t) power spectrum of the collision for the spinless reaction partners in 
the entrance and exit channels. Unlike the cross section energy autocorrelation function in 
the previous Section, the time power spectrum will be studied for a fine angular resolution. 
The time power spectrum is given by the Fourier component of the amplitude energy au- 



tocorrelation function \2l\\ . For t <C Ti/D, i.e. for Y/D 1, the continuous spectrum 



n 



approximation is valid and we have [1(3, [22] 



P{t,6) oc iJ(t)exp(-rt/^)^[iy(J)iy(J0] 1/2 exp[z($-^)(J-J0-/5|^-^1^]^jW^J'W- 

jj' 

(6) 

Here, Pj(0) are Legendre polynomials, and the Heaviside step function H(t) signifies that 
the intermediate complex cannot decay before it is formed at t = 0. 

In Fig. 2 we present P(t, 6) for three moments of time and for the two different sets of 
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the parameters for which the (7(e) 's were calculated in the previous Section (Fig. 1). The 

first set is: T = 0.15 MeV, htu = 0.75 MeV, I = 36, $ = 0, 3 = 0.03 MeV, d = 1. For the 

second set we have different values of 3 — 0.1 MeV, d = 5 while the rest of the parameters 

□ 

is unchanged. For the reason discussed in 10J the time power spectra in Fig. 2 are scaled 
with the Pdiag(t, 0), which is given by Eq. (6), where only the spin diagonal terms J = J' are 
taken into account. Such a spin diagonal approximation corresponds to the limit of quantum 
chaotic scattering and random matrix theory |l|]. Accordingly, deviation of the scaled time 
power spectra in Fig. 2 from a constant unity is a quantitative measure of the deviation of 
the collision process from the universal limit of the quantum chaotic scattering theory [l|. 

Fig. 2 illustrates a rotation of the two wave packets towards each other. As the wave 
packets rotate they also spread - the bigger 3 the faster the spreading. One observes that, 
for 3 = 0.03 MeV and d — 1, the contrast of the interference fringes, due to the interfer- 
ence between the near-side and far-side amplitudes 10(, is very strong. These interference 
fringes is a manifestation of Schrodinger cat states in highly excited quantum many-body 

1 L 

systems [22j|. On the contrary, for 3 = 0.1 MeV and d = 5, the contrast of the inter- 
ference fringes is greatly reduced indicating a quantum-classical transition in the collision 
process flo| . 

Our approach shows that the complicated many-body collision problem can be accurately 
represented by the simple picture of a weakly damped (3 « Tioj) quantum rotator. This 
picture was obtained without introducing any collective degrees of freedom of the interme- 
diate complex, such as its deformation and spatial orientation. The introduction of those 
degrees of freedom is known to be a successful approximation for very low, closed to 
Yrast line, intrinsic excitations of the intermediate complex. Yet, in our case of high in- 
trinsic excitations (> 15 MeV), the collective degrees of freedom acquire large spreading 



widths 



241 ]. r spr >> 3, r, and by consequence they decay much faster than the average 



life-time of the intermediate complex. 

Notice that P(t, 6) can be obtained from the data for excitation function fluctuations for 
binary collisions, for fine energy and angular resolutions, provided the relative contribution 
of direct processes is significant (> 70%) [25j. The latter is usually the case for heavy-ion 
elastic and inelastic scattering. Experimentally, fine energy 25J and angular [26j resolutions 
required for the determination of P(t, 6) are routinely achievable for heavy ion collisions j^ . 
Therefore, a reliable determination of the phase relaxation width 3, which is an important 
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new energy scale in quantum many-body systems 




is experimentally possible. 



IV. CONCLUSION 

We have discussed the effects of slow phase relaxation, spin off-diagonal .S-matrix correla- 
tions on the cross section energy oscillations and the time evolution of the highly excited 
intermediate system formed in the 24 Mg+ 28 Si collision. Such quasiperiodic energy oscilla- 
tions were observed experimentally. The effects of coherent rotation of the hyperdeformed 
~ (3 : 1) intermediate complex, Schrodinger cat states and quantum-classical transition 
have been revealed for the 24 Mg+ 28 Si heavy ion scattering. 

L.B. and S.Yu.K. acknowledge with gratitude financial support from the projects IN- 
101603 (DGAPA-UNAM) and 43375 (CONACyT). 
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FIG. 1: Experimental (dots) and calculated C(e)'s for 24 Mg+ 28 Si elastic and inelastic scattering. 
Dashed lines are obtained with d = 1 and (3 = 0.03 MeV, and dashed-dotted lines with d = 5 and 
P = 0.1 MeV (see text). Dotted lines are Lorentzians with T = 0.085 MeV. 
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FIG. 2: The time power spectra for the 24 Mg+ 28 Si scattering obtained for the three moments of 
time with T being a period of one complete revolution of the intermediate complex. Solid thin 
lines are obtained with d = 1 and (3 = 0.03 MeV, and solid thick lines with d = 5 and (5 = 0.1 MeV 
(see text). 



